Determination of spin Hamiltonian in the Ni4 magnetic molecule 
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Magnetic excitations in a Ni4 magnetic molecule were investigated by inelastic neutron scattering 
and bulk susceptibility (xbuik) techniques. The magnetic excitation spectrum obtained from the 
inelastic neutron scattering experiments exhibits three modes at energy transfers of hw = 0.5, 
1.35, and 1.6 meV. We show that the energy, momentum, and temperature dependences of the 
inelastic neutron scattering data and Xbuik can be well reproduced by an effective spin Hamiltonian 
consisted of intra-molecule exchange interactions, a single-ionic anisotropy, biquadratic interactions, 
and Zeeman term. Under a hydrostatic pressure, the bulk magnetization decreases with increasing 
pressure, which along with the biquadratic term indicates spin-lattice coupling present in this system. 

PACS numbers: 75.50.Xx 



I. INTRODUCTION 

Magnetic moleculesii^ provide us an excellent opportu- 
nity to study quantum behaviors in magnetism, because 
their effective spin Hamiltonians can be exactly diago- 
nalized and compared to experimental data. Previous 
studies on magnetic molecules such as Mni2 and Feg re- 
vealed superparamagnetic behaviors and tunneling be- 
tween quantum mechanical statesi^i^ These phenomena 
were explained by treating these molecules as isolated en- 
tities with strong single-axis anisotropy^ and transverse 
fourth-order term.^ It was also found that the V15 clus- 
ter exhibits a hysteresis loop with dissipative spin rever- 
sal in pulsed field magnetization measurements, which 
was explained by Landau-Zener transition and phonon- 
bottleneck effect.^ 

Ni4 is another intriguing molecular magnet because the 
antiferromagnetic Ni^+ (s = 1) ions form a tetrahedron,^ 
which may lead to geometrically frustration. The full 
chemical formula of the Ni4 cluster is [Moi203o(^2- 
OH)ioH2{Ni(H20)3}4]-14H20, and its crystal structure 
is shown in Fig. [TJ four Ni^+ ions form a slightly dis- 
torted tetrahedron, and tetrahedra directing oppositely 
in the c-axis are arranged alternately. Distances between 
Ni^+ ions within a cluster are 6.69, 6.70, 6.62, and 6.60 A, 
whereas the shortest distance between Ni^+ ions that be- 
long to different clusters is 7.15 A. Previous bulk prop- 
erty measurements using bulk susceptibility, high field 
magnetization, electron paramagnetic resonance, optical 
conductivity, and magneto-optical response^"— showed 
that the dominant interaction in the system is antiferro- 
magnetic, which is due to the superexchange interaction 
through the Ni-0-Mo-O-Ni bonds.? The most interesting 
property of the Ni4 nanaomagnet is an adiabatic change 
with non-equidistant steps observed in the magnetiza- 
tion measurements as a function of external magnetic 
fieldi2r2i This was explained by a model Hamiltonian 
that consists of a field-dependent exchange parameters. 



a single- ion anisotropy, and a biquadratic interaction.— 
The spin Hamiltonian should be determined by a more 
direct tool than the bulk property measurements to check 
its validity. 

We have performed inelastic neutron scattering mea- 
surements on a powder sample of the Ni4 molecule to 
investigate the energy {fvjj), momentum (Q), and tem- 
perature (T) dependences of the magnetic excitations. 
Discrete excited levels were observed at 0.5, 1.35, and 
1.6 meV, and Q-dependences of each mode have peaks at 
0.6 and 1.6 A"-'^. By analyzing the inelastic neutron scat- 
tering data, we determine the spin effective Hamiltonian 
of the spin Ni4 cluster without an external field; a model 
Hamiltonian consists of intra-molecule exchange interac- 
tion, single-ion anisotropy, and biquadratic interaction.— 
The existence of the biquadratic interactions suggests 
strong spin-lattice coupling in the system, which is con- 
sistent with the suppression of the bulk magnetism by an 
application of a hydrostatic pressure. 




FIG. 1: (Color online) Schematic view of the crystal struc- 
ture of the Ni4 magnetic molecule. Solid circles represent 
Ni^"*" ions. Solid and dashed lines represent the exchange and 
biquadratic interactions, and their parameters, J, J', j, and 
J , are defined in Eq. ((2|. 



II. EXPERIMENTAL DETAILS 

A 4.5 g deuterated powder sample was prepared using 
the procedure described in Ref. |8|. Our prompt-gamma 
neutron activation analysis showed that about 50% hy- 
drogen was substituted by deuterium. A small amount 
of the sample of 200 mg was used for bulk magnetization 
using a SQUID magnetometer at the ambient pressure as 
well as under the hydrostatic pressure up to 0.92 GPa. 
For the pressure experiment, 20 mg sample was put into a 
teflon cell which was then filled with Daphne oil and was 
set to a piston-cylinder devicej^ A reference sample of Sn 
was also put in the cell, and the transition temperature 
of Sn was used to determine the hydrostatic pressureji2, 
Background for the pressure experiments was measured 
using the empty pressure device, and subtracted from the 
data. 

The remaining 4.3 g sample was used for two sets 
of neutron scattering experiments. The first set of the 
measurements was performed on the cold-neutron triple- 
axis spectrometer SPINS at the NIST Center for Neu- 
tron Research. A vertically focusing pyrolytic graphite 
(PG) monochromator and a horizontal focusing PG an- 
alyzer were used to increase the sensitivity of the mea- 
surements. Energy of the scattered neutrons was fixed to 
be Ef — 3.0 meV, resulting in an instrumental resolution 
of 117 fj.eV (FWHM, or full width at half maximum) at 
the elastic position. Energy of the incident neutrons was 
changed to measure the scattering intensity as a func- 
tion of energy transfer, huj. The energy resolution at 
fkjj = 0.5, 1.35, or 1.6 meV is estimated to be 141, 189, or 
212 iieV (FWHM), respectively^^ Higher order contam- 
inations were eliminated using a cooled Be filter placed 
after the sample. The second set of experiments was per- 
formed at the cold neutron triple-axis spectrometer HER 
at the JRR-3M research reactor with Ef = 5.0 meV. A 
vertically focusing monochromator and a double-focusing 
(i.e., both horizontal and vertical focusing) analyzer were 
employed to increase the sensitivity. The nonmagnetic 
background was measured at 30 K and subtracted from 
the low temperature data to obtain the magnetic scat- 
tering intensity, I{Q,hiu), 

/(Q,M = 5]^^(<5„,^-Q"g^/g2)F2(Q) 

a,/3 a,b i,f 

X p, < zi^^-'Q-"-"!/ X /|5,V^-^1» >(i) 

X S{Ei - Ef + huj) 

where \i > (|/ >) is the initial (final) eigenstate. 
Pi is the Boltzmann factor for the state \i > {pi = 
n-ie"^'/'"'^-^/^ rij-e"^^/'^^^ where rii represents degen- 
eracy of \i >), Ei (Ef) is the initial (final) energy of the 
system. So, 6 and Va^t are the spin operator and position 
of the Ni^+ ions at site a, b in the molecule, respectively, 
a and /? represents the vector component of S, and F{Q) 
is the magnetic form factor of the Ni^+ ion J^ 
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FIG. 2: (Color online) (a) The fiij-dependences, I{hjjS), mea- 
sured at Q = 0.5 A"^ and T = 2.0 and 30 K. The fej- 
dependences of the magnetic intensities at (b) 2.0 K ob- 
tained by I{huj, 2.GK) - 7(fe,30K) and at (c) 4.0 K by 
/(ftw, 4. OK) —I{huj, 30K), respectively. Solid and dashed lines 
are described in the main text. 



III. RESULTS AND DISCUSSION 

Figure ^a) shows the ^loj-dependence of the neutron 
scattering intensities measured at Q = 0.5 A~^ at two 
different temperatures of T = 2.0 and 30 K. The strong 
sharp peak centered a,t hoj — meV is mainly due to the 
incoherent scattering from hydrogen in the sample. At 
T = 2.0 K, in addition to the strong incoherent scatter- 
ing, three distinct inelastic excitations exist centered at 
around hoji = 0.5, huj2 = 1.35, and htos = 1.6 meV. 
At T = 30 K, the three inelastic peaks broaden and 
become indistinguishable with the incoherent scattering. 
We take the T = 30 K data as the background and sub- 
tract it from the T = 2 K data to obtain the magnetic 
energy spectrum. The result, /(2K) — /(30K), is plot- 
ted in Fig.^b). In order to investigate the temperature 
dependence, we have performed the same measurements 
at 4.0 K. /(4K) - /(30K) measured at Q = 0.5 A'^ is 
shown in Fig.I^Jc). The intensities of the three peaks be- 
come weaker at 4.0 K, indicating that these three peaks 
are magnetic. In order to study T-dependence further, 
we have measured T-dependence of the huji mode and 
bulk susceptibility (Xbaik)- As Fig. [3] shows, upon cool- 
ing the hui mode slowly appears at low temperatures 
and rapidly increases below 5 K, which coincides with 
the downturn in Xbaik- This supports that the excita- 
tions observed in the neutron scattering spectra are due 
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FIG. 3: (Color online) (a) T-dependence of the hwi — 
0.5 meV mode measured at Q = 0.55 A~^. Solid line is 
the Boltzmann factor for the ground state, po explained in 
the main text. Dashed line represents background, (b) Ex- 
perimental (dots) and model calculated (line) T-dependence 
of the bulk susceptibility, Xbuik, under an external magnetic 
field of B = 0.01 T. 
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FIG. 4: (Color online) Q-dependences of (a) elastic, (b) huji — 
0.5 meV, (c) hiU2 = 1.35 meV, and (d) huu^ — 1.6 meV modes. 
Solid lines are explained in the text. 



to a development of antiferromagnetic correlations. 

We have also measured the Q-dependences of the three 
excitations at T = 2.0 or 0.7 K. As shown in Figs.lUJb)- 
IDJd), the 0.5, 1.35, and 1.6 meV excitations are peaked 
at Qo — 0.6 and 1.6 A"-'^, while no clear magnetic in- 
tensity is observed in the Q-dependence at the elastic 
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FIG. 5: (a) Energy levels and eigenstates of the Ni4 magnetic 
molecule for "Ho ~ — J^^Si • Sj. The subscript, a ~ e, rep- 
resents the degeneracy of each state. Some eigenstates are 
described in Table HI (b) Some of the energy levels and eigen- 
states of the Ni4 spin cluster with H written in Eqs. ([2]) and 
Q. Some eigenstates \n > are described in Table [III 



channel [Fig. UJ^a)]. All three excitations weaken at high 
Q, confirming that they are magnetic in origin. The Q- 
dependence of the three excitations are similar to each 
other, indicating that they have the same origin. The 
Q-dependence of the intensity has the information about 
what is the magnetic entity in real space, which will be 
discussed later in detail. 

Following the previous study,~ we assume that the 
spin Hamiltonian of the Ni4 spin cluster consists of the 
exchange interaction, the single ion anisotropy, the bi- 
quadratic interaction, and the Zeeman term; 



'H= - Y. J^JS^ ■ ^3 + Dj2{e^ ■ S,)^ (2) 

4 4 

i=ij=l i=l 

where Jy, D, and jtj are the parameters of the ex- 
change interaction, the single ion anisotropy, and the bi- 
quadratic interaction, respectively. By considering the 
crystallographic symmetry in the Ni4 moleculcfi^ there 
are two different values for Jij and jij as described by 
the sohd and dashed lines in Fig. (TJ and we define J, J', 
j and j' using Jij and jij as J = J12 = J13 = J14, 
J' = ^3 = -^24 = -^34' 3 = 312 = ji3 = ji4, and 
/ = J23 = i24 = ik- 9 and ^b represent the geo- 
metric factor and the Bohr magneton, e^ describes a 
local anisotropic axis; considering the geometrical frus- 
tration, e^ points radially outward from the center of the 
tetrahedron.^" 

If the Ni4 cluster has uniform exchange couplings only, 
the spin Hamiltonian will be simply written as T-Lq = 
— J^Si ■ Sj. The ground state of T-Lq is a triply de- 
generate state with zero total spin, l^totai = >i^3. 
All the eigenstates of "Ho can be easily calculated, some 



TABLE I: Some of low energy eigenstates of the Ni4 magnetic molecule assuming it has the simple hamiltonian Ho described 
in Fig. [S{a). The eigenstates are described by the linear combinations of the ^-components of the four Ni^+ [s = 1) spins 
\s1S2S3S4 >. t) 0, and 4. represent s^ — 1, 0, and -1, respectively. For simplicity, the terms whose coefficients are less than 0.2 
were not written here for the linear combinations, but they were included in our calculations. 
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5total = 2 >3 



-0.25 



0.23| ttii> -0.25100 t4-> -0.25|00 4,t> +0.24| t 00 t> -0.22|0 t t> +0.23| Utt> -0.29| t4-tt> 
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-0.22J0 t ;> +0.22I tt i> +0.25I U 00 > +0.25| tt 00 > -0.27|00 t > -0.22|0 t t> 
+0.26] ti t> 



2>2 



-0.22|00 tt> +0.27I ttt4-> +0.28I ttt > "0.22| t tt> +0.23| t4.tt> 



Stotal — 2 >2 



0.23|00 tt> -0.27|0 ttt> +0.21|32 > +0.22| tt 00 > +0.22| tt 00 > -0.27| t t > 
-0.21| ttOt> +0.21| tttt> 



2>2 



0.23|00 tt> +0.24| t tt> -0.23| ttU> +0.35| ttt4-> -0.23| ttt4-> +0.2| tt4-t> 



Ototal — 2 >2 



0.4| UU> -0.26| t 00 t> -0.32|0 tt > +0.24| t t > +0.27| tt 00 > 



2>2 



0.24| t 00 t> -0.27| ttt^> +0.2910 tt4-> +0.23| t t4-> 
0.22| ttit> +0.23] t tt> -0.35] ttit> 



-0.26| tOtO> -0.23|0t00> 



>5total — 2 >l 



-0.22| t tt> -0.38| ttii> +0.23J t 00 t> +0.32| t t > -0.37| tt 00 > 



Stotal — 2 >l 



0.27| t 00 t> -0.21|0 iU> -0.21| t U> 
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-0.28| t t > -0.2|0 U0> +0.21|0 t 00 > 



-0.310 tti> +0.2| t 00 t> +0.23I t U> +0.22| ttt~l> +0.35| tti > 

+0.25I ttit> +0.2| t 00 t> 
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5total = 1 >17 



1>1 
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~0.22j ttit> -0.24J t tt> +0.331 ttit> +0.3| tt t> 



-0.28] tttO > -0.21| tOtO> 



l>i 



Stotal — 1 >1 



-0.31| t ti> +0.23I tt t> +0.21|16 > -0.21| tttt> 
+0.43I ttt0>+0.21| ttit> 



-0.27| t ti> -0.25|0 t 00 > 



Stotal — 1 >1 
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-0.22!0tt0> 



•S'total — >2 



-0.311 tttt> +0.33100 tt> -0.351 tttt> +0.331 tt 00 > -0.3310000 > +0.331 tt> -0.351 tttt> 
+0.33100 tt> -0.31| tttt> 



0>i 



0.291 t 00 t> -0.2910 t t> - 
-0.291 t t > +0.27I tttt> 



0.31| tttt> +0.27I tttt> -0.291 t t > +0.29|0 tt > +0.29|0 tt > 
-0.31| tttt> -0.29| t 00 t> +0.291 t 00 t> 



of which at low energies are illustrated in Fig. [5{a) and 
listed in Table HI For the general H with non-uniform 
exchange J^, single-anisotropy D, and biquadratic jij, 
the eigenstates can be written as linear combinations of 
the eigenstates of Hq. For instance, the ground state of 



H, denoted by |0 >, becomes a singlet state that is a 
linear combination of the triply degenerate ground state 
I -S'total = >i^3 of Hq. The other two linear combina- 
tions of [S'total = >i.,.3 gain energy for H. These states 
and other low-energy excited states, and their energies 



TABLE II: Some of low-energy eigenstates and energies of the Ni4 magnetic molecule for the Hamiltonian Ti of Eqs. ^ and 
(|3]). The eigenstates of H can be described by hnear combinations of the states of Ho shown in Table H) For simpUcity, the 
terms whose coefficients have an absolute value less than 0.3 were not written here for the linear combinations, but they were 
included in our calculations. 



£„ - Eo (meV) 



n > 



1.930 



24 >~ (-0.76 + 0.02i)|5'totai = 1 >i7 +(-0.4 - 0.03i)|gtotai = 1 >i 



1.930 



23 >~ (-0.08 + 0.4i)jS'totai = 1 >i7 +(0.22 - 0.73i)|5'totai = 1 >i 



1.869 



22 >:::; (0.65i)|5'totai = 2 >i7 +(-0.02 - 0.3i)|gtotai = 2 >29 



1.869 



21 >~ (-0.01 + 0.33J)|S'totai = 2 >i8 +(-0.02 - 0.35J)|S'totai = 2 >2i +(0.01 + 0.33i)|Stotai = 2 >2 

+ (0.02 - 0.5l)|Stotal = 2 >27 +(0.03 - 0.32j)|5'total = 2 >28 



1.862 



20>~(0.53-0.23J)|Stotai = 
+(0.34-0.16i)|gtotai = 



■ 2 >19 +(-0.31 + 0.14i)|S'total = 2 >21 +(-0.32 + 0.14i)|Stotal = 2 >23 
2>28 



1.814 



19 >~ (-0.49 - 0.01J)|5'totai = 2 >2i +(0.63 + 0.02J)|S'totai = 2 >23 +(-0.31 + 0.05i)|Stotai = 2 >24 



1.814 



18>~(0.36-0.08i)|5totai = 
+(0.36-0.07i)|gtotai = 



2 >17 +(0.3 - 0.1J)|S'total = 2 >19 +(0.31 - 0.11i)|5total = 2 >20 
2>24+(-0.49 + 0.13i)|gtotal = 2 >30 



= 1 >i +(0.2 + 0.31i)|5totai = 1 >3 +(0.04 - 0.37J)]5'totai 
: 1 >5+(0.34 + 0.34^)|5'total = 1 >6 



1.796 



17>~(-0.46 + 0.15J)|S'totai 

+ (-0.3 + 0.31i)|Stotal = 



1 >4 



16 >- (-0.21 + 0.3H)|5'totai = 1 >i +(0.5 + 0.39i)|gtotai = 1 >i +(-0.23 + 0.32i)|5'totai = 1 >3 

+ (0.37 - 0.24i)jgtotai = 1 >3 +(0.3 - 0.08^)|gtotal = 1 >4 
r, +("-0 5.^ - 31j,';I,S'*„*„, = 1 ^« 



1.726 



1.726 



15>~(0.19 + 0.57i)|Stotai = l>2 



1.669 



14 >~ (-0.32 + 0.56i)|5'totai = 1 >5 +(-0.55 - 0.31z)|gtotai = 1 >6 



1.669 



13 >~ (-0.49 + 0.34J)|S'totai = 1 >i +(-0.04 - 0.42z)|Stotai = 1 >3 +(0.42 + 0.34i)|5'totai = 1 >4 



1.632 



12>~(0.21+0.35i)|Stotai = 
+(-0.03-0.49i)|gtotai 



1 >3 +(0.44 - 0.26i)j5'totai = 1 >4 +(0.44 - 0.02i)|5'totai = 1 >5 

= 1>6 



1.424 



11 >~(-0.16-0.48J)|S'totai 
+(0.12 + 0.37i)|Stotai = 



= 1 >8 +(-0.47 + 0.2i)|S'total = 1 >9 +(-0.3 + 0.04i)|S'total = 1 >io 
1 >14 



1.383 



10>~(-0.4-0.09i)|5'totai 

+(-0.37-0.12i)|gtotai 



= 1 >7 +(-0.16 + 0.51i)|5total = 1 >8 +(-0.29 + 0.24i)|5'total = 1 >10 
= 1 >13 



1.383 



9>~ (0.33-0.47i)|5'totai = 
+(0.34 + 0.18i)|Stotai = 



1 >9 +(0.34 + 0.1i)|5total = 1 >10 +(0.32 + 0.18i)|Stotal = 1 >14 
1>15 



1.306 



>~(0.37-0.15i)|5'totai = 
+(-0.2 + 0.25i)|gtotai 



1 >7 +(-0.42 + 0.03i)|Stotai = 1 >9 +(0.32 + 0.33J)|S'totai = 1 >io 
l>ii 



1 >i2 +(-0.05 + 0.52j)|5'totai = 1 >i3 +(0.25 + 0.28i)|Stotai 

= 1>15 



1.243 



7>~(0.53 + 0.07i)|5'totai = 
+(-0.29 + 0.18i)|Stotai 



1 >14 



1.243 



6 >~ (0.31 - 0.43^)|5'total = 1 >i2 +(-0.47 - 0.23^)|5'total = 1 >i3 +(-0.31 + 0.17^)|gtotal = 1 >i5 



0.714 



5 >~ (0.55 - 0.06i)|5'totai = >i +(-0.49 + 0.07i)|5totai = >2 +(-0.02 - 0.65i)|S'totai = >3 



0.714 



4>~ (0.36-0.33z)|5't, 



:0>i +(-0.33 + 0.29i)|gtotai =0>2 +(0.44 + 0.6i)|5'totai = >3 



0.531 



3 >~ (-0.6H)|5'total = 1 >7 +(-0.34i)lgtatal = 1 >8 +(0.02 + 0.57i)|gtotal = 1 >1 



0.531 



2 >~ (0.78 + 0.26i)|5't, 



l>i 



0.512 



1 >~ (0.28 + 0.5z)|5'totai = 1 >i2 +(-0.26-0.47i)|gtotai = 1 >i4 +(0.27 + 0.49i)|gt, 



= 1 >15 







>~ (-0.33 - 0.58^)|5'total = >i +(-0.37 - 0.65i)|gtotai = >2 



were obtained by the exact diagonalization of 'K for sev- 
eral different sets of values for the parameters J, J', Z?, 
j, and i' . The optimum parameters were determined by 
comparing both the calculated energies and intensities of 
the allowed transitions between the states with the ob- 
served energies and intensities of the excitation modes, 
hw\, fiW2, and }vjJ^\ 



(3) 



This result clearly show the biquadratic interactions 



Jlk^ = 


-3.69(3) K 


J'/k^ = 


-3.19(2) K 


D/kB = 


-2.47(2) K 


3/kB = 


-0.11(1) K 


j'/kB = 


1.52(1) K. 



eigenstates of H for the optimal parameters are listed 
in Table. HT] The fiw-dependence of neutron scattering 
intensities for the three low-energy excitations were cal- 
culated by Eq. ([1]), averaged for powder, and convoluted 
with Gaussians. The fitting results shown as solid lines 
in Figs. [2{b) and ^c) reproduce our data well. The Q- 
dependences of the excitations can be also reproduced 
by the model as shown in Figs. |lUb)-[5l^d). On the other 
hand, the calculated /iw-dependences using the previ- 
ous parameters reported in Ref. [ifl (J/fce = —3.2 K, 
J'/kB = -3.1 K, D/kB = -1-0 K, j/kB ^ 1.6 K, and 
j'/^B = K) cannot fit the experimental ftw-dependences 
at all [see the dashed lines in Figs.IUJb) and[2jc)]. 

T-dependence of the hioi mode can be reproduced very 
well by the Boltzmann factor for the ground state, po — 



present in the Ni4 cluster. The resulting low-energy ^/^j^jC ^j/*^^'^, as described by solid hne in Fig. [HJa) 
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FIG. 6: (Color online) T-dependence of the bulk magnetiza- 
tion, Mbuik, measured with B = 0.1 T under a hydrostatic 
pressure of P = 0, 0.38, 0.54, and 0.92 GPa. P-dependences 
of (b) the majcimum value of Mbuik {M^^) and of (c) the 
peak temperature (T*). 



T-dependence of Xbuik can also be reproduced well by our 
model as shown in Fig. ^h). In the calculation of Xbaik, 
the geometrical factor was taken to be 5 = 2.22.—"— The 
background due to isolated Ni^+ ions was estimated by 
the Curie law and added to the calculated Xbuik- 

The existence of the biquadratic terms in H indicates a 
spin-lattice coupling in this system. ^^ In order to confirm 
this, we have measured the bulk magnetization, Mbuik, 
under a hydrostatic pressure (P). Figure [Slja) shows 



A^buik obtained with B = 0.1 T under a hydrostatic 
pressure of P = 0, 0.38, 0.54, and 0.92 GPa. As the 
pressure increased, the maximum value of Mbuik (-^bufk) 
decreases [Fig. H^b)], while the peak temperature (T*) 
increases [Fig. |6l[c)]. These results represent that P en- 
hances antiferromagnetic correlations by shortening the 
distance between the Ni^"*" ions, which may explain the 
existence of the biquadratic term in this system. 



IV. CONCLUSION 

We have determined the effective spin Hamiltonian, 
T-L, in the deuterated Ni4 magnetic molecule by using in- 
elastic neutron scattering and exact diagonalization tech- 
niques, fiw-, Q-, and T-dependences of neutron scat- 
tering intensities due to the low-huj excitations centered 
at 0.5, 1.35, and 1.6 meV as well as T-dependence of 
the bulk susceptibility can be well accounted by Ti con- 
sisted of exchange interaction, single-ion anisotropy, bi- 
quadratic interaction, and Zeeman term. 
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